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The galaxy M82 has long been considered a promising target for very-high-energy (VHE) gamma-ray observa- 
tions because of the compact starburst region in its core. Theoretical predictions have suggested it should be 
detectable by ground-based imaging Cherenkov telescopes like VERITAS and that a detection would have im- 
plications for the understanding of the origin of cosmic rays. M82 was observed with the VERITAS array during 
the 2007-2009 observing seasons. With an exposure of 137 hours, VERITAS was able to detect a gamma-ray 
signal at the 5<r level. This marks the discovery of gamma rays not only from M82 but also from the new source 
class of starburst galaxies. The observed flux from M82 is (3.7 ± 0.8 a tat ± 0.7 sys t) x 10~ 13 photons cm -2 s~ 1 
above an energy threshold of 700 GeV, which corresponds to 0.9% of the Crab Nebula flux. The differential 
energy spectrum is a power law with a photon index T = 2.5 ± 0.6 s tat ± 0.2 sys t. Both the flux and the photon 
index are close to recent theoretical predictions. The VERITAS data indicate a strong correlation between the 
star-formation activity and the cosmic-ray production in M82. 



1. Introduction 

The existence of cosmic rays near the Earth has 
been well established ever since their discovery in the 
early 20th century but their origin has eluded us 
ever since. Although there are no direct measurements 
of cosmic rays beyond Earth, there are indirect ob- 
servations which give strong evidence for cosmic rays 
permeating the entire Galaxy and most significantly 
in the Galactic plane. This comes from the diffuse 
gamma-ray emission seen by the Fermi Large Area 
Telescope (Fermi-LAT), interpreted as mainly com- 
ing from cosmic-ray ions interacting with interstellar 
gas producing neutral pions which subsequently decay 
into gamma rays, although there is also a contribution 
from cosmic-ray electrons upscattering ambient pho- 
tons to gamma-ray energies. 

The leading hypothesis is that winds and super- 
novae of massive stars provide the major accelera- 
tion sites of Galactic cosmic rays, and observations 
of supernova remnants, such as the shell-type super- 
nova remnant RX J1713. 7-3946 [2|, are suggestive of 
this scenario. The detection of the Large Magellanic 
Cloud (LMC) with the Fermi-LAT Q shows gamma- 
ray emission coincident with the star-forming region 
30 Doradus in the LMC. If massive star winds and 
supernova remnants accelerate cosmic rays then star- 
forming regions like this should indeed be emitting 
gamma rays. 

The bright galaxy M82 is located at a distance 
of about 3.4 Mpc in the direction of the constella- 
tion Ursa Major Ogoing interactions with nearby 
galaxies, including the spiral galaxy M81 (see [j|), 
have formed a compact active starburst region, about 
1000 light yrs across, in the center of M82. In this re- 
gion, stars are being formed at a rate approximately 
ten times faster than in an entire "normal" galaxy like 



the Milky Way and the rate of supernovae is 0.1-0.3 
yr" 1 [1,0]. Observations of the central region of M82 
at radio frequencies, interpreted as synchrotron radi- 
ation of cosmic-ray electrons spiraling in the galactic 
magnetic fields, suggest a very high cosmic-ray energy 
density, about two orders of magnitude higher than in 
the Milky Way @]. The mean molecular gas density 
is also high, about 150 particles per cubic centimeter 
or about 10 9 solar masses total in the starburst region 
0- 

Because of the aforementioned properties of the 
starburst core, M82 has long been considered a proba- 
ble gamma-ray source. Neither the Energetic Gamma- 
Ray Experiment Telescope (EGRET) [Kj onboard 
NASA's Compton Gamma-Ray Observatory nor pre- 
vious ground-based gamma-ray observatories, such as 
the Whipple 10m telescope \V\\ . detected gamma-ray 
emission from M82. Upper limits on the gamma-ray 
flux were set by EGRET at 4.4 photons cm~ 2 s _1 for 
E > 100 MeV and by Whipple at about 10% of the 
flux from the Crab Nebula for E > 100 GeV. The lat- 
ter limit is above the sensitivity of the Very Energetic 
Radiation Imaging Telescope Array System (VERI- 
TAS). 



2. Observations and Analysis 

VERITAS fl2f is an array of four imaging atmo- 
spheric Cherenkov telescopes located at the basecamp 
of the Fred Lawrence Whipple Observatory near Tuc- 
son, Arizona. The array has been fully operational 
since mid-2007 and has the sensitivity to detect a 
point source with 1% of the steady Crab nebula flux 
in less than 50 hours 1 . Observations cover the energy 



1 The integral flux sensitivity above 300 GeV was improved 
by about 30% with the move of one telescope in the summer of 
*see http://veritas.sao.arizona.edu/conferences/authors7icrc2009 2009, which corresponds to a decrease in the time required to 
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Figure 1: Sky map showing the measured excess of VHE 
gamma-ray-like events above the estimated background. 
The excess is centered on the location of M82. Each pixel 
in the map contains the excess in a circular region of 
radius 0.1° and the background in the pixel is estimated 
using an annulus centered on the pixel. The spatial 
distribution of the excess is consistent with a point-like 
source. The circle is a representation of the VERITAS 
point spread function (68% containment) for individual 
gamma rays. The black star denotes location of the core 
of M82. (Figure from [jj].) 

range from 100 GeV to beyond 30 TeV with an energy 
resolution of 15-20% above 300 GeV, and an angular 
resolution per gamma-ray photon of 0.1° at 1 TeV. 

The observations of M82 were made in the period 
between January 2008 and April 2009. The total ex- 
posure amounts to about 137 hours of quality-selected 
live time, i.e. time periods of astronomical darkness 
and clear sky conditions, which is the deepest expo- 
sure taken with VERITAS to date. By observing the 
source offset from the center of the field-of-view, si- 
multaneous estimation of the background was made 
possible [l3|. The source was observed at a mean 
zenith angle of 39° . 

The data analysis was performed using the standard 
VERITAS analysis tools Q 

using event-selection cri- 
teria that were optimized a priori for a low-flux hard- 
spectrum source. A total excess of 91 gamma-ray- 
like events (see the sky map in Figure [T|) were de- 
tected above the estimated background (267 back- 
ground events). The excess corresponds to a post- 
trials significance of 4.8cx. The observed gamma-ray 
flux above the 700 GeV energy threshold of the anal- 



detect a point source with 1% of the steady Crab flux to about 
30 hours. This did not have any impact on the M82 data set. 
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Figure 2: Time-averaged gamma-ray flux during each 
monthly observation period. The integral flux was 
calculated with the assumption of a power-law spectrum 
with a photon index V = 2.5. The error bars represent la 
statistical errors only. The horizontal bars reflect the 
actual range of dates for each exposure. The dashed line 
shows the fit of a constant value, indicating a steady 
steady gamma-ray flux. (Figure from supplementary 
material to [l6T|.) 

ysis is (3.7 ± 0.8 sta t ± 0.7 syst ) x 10~ 13 photons cm~ 2 
s _1 and, as shown in the lightcurve in Figure [H there 
are no flux variations. From the gamma-ray flux we 
infer the gamma-ray luminosity to be 2 x 10 32 W, 
which is about 2 x 10 6 times smaller than the measured 
far- infrared (100/im) luminosity [l||. The differential 
VHE gamma-ray spectrum (shown in Figure[3]) is best 
fitted using a power-law function dN/dE oc E~ r with 
a photon index of T = 2.5 ± 0.6 s t a t ± 0.2 syst (the two 
uncertainties corresponding to statistical and system- 
atic errors respectively). 

The VHE gamma-ray flux from M82 is 0.9% that 
observed from the Crab Nebula, which makes it one 
of the weakest gamma-ray sources ever detected. Be- 
cause of the exceptionally long exposure, several tests 
were performed to ensure that systematic effects did 
not introduce a spurious signal in the data. The signal 
and the measured spectrum and flux have been veri- 
fied with independent calibration and analysis chains. 
For further details see the supplementary information 
of 0. 



3. Discussion 

VERITAS has for the first time detected VHE 
gamma-ray emission from an extragalactic source that 
is not clearly assiociated with an active galactic nu- 
cleus (AGN), in which the emission of gamma rays is 
powered by accretion onto a supermassive black hole. 
It is possible that M82 harbors a supermassive black 
hole in its center but there is at most very weak signs 
of AGN activity 

The high star formation rate in the core of M82 im- 
plies the presence of strong shocks which are expected 
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Figure 3: Differential energy spectrum (E x dN/dE, 
where E is the gamma-ray energy and N is the photon 
count) of M82 measured with VERITAS (open diamonds 
with la error bars). The VERITAS points can be fitted 
with a power-law function, dN/dE oc E~ r , where 
T = 2.5 ± 0.6 s tat ± 0.2 syst , indicated by the thick grey 
line. The flux upper limit at about 6.6 TeV is above the 
extrapolation of the fitted power law at that energy. The 
theoretical model prediction [2p| ]. of the total emission is 
given by the thin solid line. Its components are from tt° 
decays, inverse Compton scatterings and bremsstrahlung. 
(Figure from flt|.') 



to accelerate cosmic rays. These cosmic rays perme- 
ate the galaxy and produce gamma rays as they in- 
teract with interstellar gas and photon fields. Recent 
theoretical work [TH, O H(| has predicted the VHE 
gamma-ray flux from M82 based on acceleration and 
propagation models of cosmic rays in the starburst re- 
gion and the galaxy. These authors all predict fluxes 
close to that measured by VERITAS. In Figure [3] the 
model of HO] is compared with the VERITAS result. 

The cosmic-ray density in the starburst core of M82 
can be estimated using the observed VHE gamma-ray 
flux. At about 250 eV cm' 3 , the density is about 500 
times larger than the average of the Milky Way. It 
must be noted though, that the Milky Way is much 
larger in terms of volume and thus the total energy in 
cosmic rays is similar. 

Because of various energy-loss mechanisms, such as 
adiabatic cooling in the stellar winds and interactions 
with interstellar gas, the lifetime of cosmic rays in M82 
is of the order of 1 Myr. This is lower than the typical 
lifetime of cosmic rays in the Milky Way by about a 
factor of 30. In order to maintain the same level of 
energy content in cosmic rays, the power input into 
cosmic rays must be correspondingly higher. Interest- 
ingly, the supernova rate is about 30 times higher in 



the starburst region of M82 than in the entire Wilky 
Way. The VERITAS detection of VHE gamma rays 
from M82 shows a correlation between cosmic-ray ac- 
celeration and massive-star formation and provides an 
important piece in the effort to reveal the origin of 
cosmic rays. 

The observed VHE gamma-ray emission includes 
contributions from both hadronic (protons and heav- 
ier ions) and leptonic (electrons and positrons) in- 
teraction channels. The spectra from these channels 
are quite different in the VHE regime. Cosmic-ray 
ions produce gamma rays in collisions with interstellar 
gas. In these collisions, unstable secondary particles, 
mostly pions, are produced which subsequently decay. 
The electrically neutral pions decay directly into two 
photons and the charged pions decay into neutrinos 
and muons, which then decay into more neutrinos and 
electrons or positrons. The latter can then in the pres- 
ence of a magnetic field emit synchrotron radiation in 
the radio or infrared bands. 

If one assumes that the electron population is 
strictly secondary electrons from the above mentioned 
process, then an observation of the synchrotron radio 
emission sets an upper limit on the gamma-ray flux 
from cosmic-ray ions. The observations of M82 at 32 
GHz frequency [2lJ set the limit at 2.5 x 10~ 9 photons 
cm~ 2 s _1 at a photon energy of 20 GeV (assuming the 
magnetic field is not much weaker than the current 
estimate of 8 nT). This energy is below the detec- 
tion capabilities of VERITAS, but an extrapolation 
of the VERITAS spectrum (T = 2.5) to this energy 
would exceed the limit by about a factor of two. This 
means the spectrum at TeV energies has a slightly 
lower photon index than the VERITAS best-fit spec- 
trum or the gamma-ray emission is not predominantly 
from cosmic-ray ions. 

An alternative scenario is that the observed radio 
emission comes from primary electrons accelerated in 
the starburst region. These electrons would also in- 
teract with ambient photon field and upscatter those 
photons to hard X-ray and soft gamma-ray energies. 
By observing the non-thermal X-ray emission it is 
then possible to constrain the electron population. 
Observations place the non-thermal X-ray luminosity 
at 5 keV photon energy, not much higher than the ob- 
served VHE gamma-ray luminosity. The X-ray data 
place a lower limit on the interstellar magnetic-field 
strength at about a third of the current estimate of 
8 nT. From this an upper limit on the absolute num- 
ber of electrons with energies around 1 GeV can be 
derived. Theoretical predictions suggest the gamma- 
ray spectrum from Compton upscattering should be a 
power law with a photon index of about 2 in the 100 
keV to 100 GeV energy range. Electrons with ener- 
gies higher than that will lose energy quickly, and the 
acceleration efficiency is decreased at a characteristic 
energy. This implies a cut-off should be introduced 
in the spectrum and identification of such could help 
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identify the underlying gamma-ray production mech- 
anism. 

Recently, the Fermi-LAT team has reported on the 
detection of high-energy gamma rays from M82 [22j ■ 
Combination of data from VERITAS and Fermi-LAT 
will further help us understand the origin of cosmic 
rays. 
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